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ABSTRACT 

The coupling between a two-dimensional, supersonic, laminar boundary layer and a flex- 
ible surface is studied using direct numerical computations of the Navier-Stokes equal, cuts 
coupled with the plate equation. The flexible surface is forced to vibrate by plane aeons c 
waves at normal incidence emanated by a sound source located on the side o ,c 
surface opposite to the boundary layer. The effect of the source excitat.on frequency on 
'surface vibration and boundary layer stability is analyzed. We find that, for frequences near 
the fifth natural frequency of the surface or lower, large disturbances are introduced ill the 
boundary layer which may alter its stability characteristics. The interaction between a sta 
two-dimensional disturbance of Tollnnen-Schlichting (TS) type with the vibrating surface s 
also studied. We find that the disturbance level is higher over the vi mating exi e sur 
than that obtained when the surface is rigid, which indicates a strong coupling “ 
and structure. However, in the absence of the sound source the disturbance level - e £ 
rigid and flexible surfaces are identical. This result is due to the high frequency of the 
disturbance which does not couple with the flexible surface. 
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Introduction 


In recent years, the demand for developing a high speed civil transport has increased. 
This has led to an increase in research activity on compressible supersonic flows, m par- 
ticular on the evolution of unsteady disturbances in a supersonic laminar boundary layer. 
One Cass of unsteady disturbances which ha, received considerable attention ■ mstab.hty 
waves in a boundary layer, i.e. eigen-modes of the compressible Orr-Sommerfeld e<iuat,ons 
obtained via linearisation around a parallel flow. When these waves are unstable sma 
disturbances cart evolve to large nonlinear disturbances as they propagate downstream m 
the boundary layer, leading to transition from laminar to turbulent flow. 

The linear stability of compressible laminar boundary layers has been stud 
sively for both subsonic and supersonic flow regimes. These studies have led to the well 
known partition of instabilities into two different classes; the viscosity dominated Cass 
known as the verity or ** mode, which is similar to the 

type waves found in low speed flows, and the acouA or HgUr modes • - . At lngh 
numbers (> 3, for adiabatic conditions), it is found that the dominant modes of inst. 
bility are the acoustic ones, and specifically the first acoustic mode known a, the «cun 
mode The linear stability of a two-dimensional flat plate boundary layer including >o i 
the vorticity and the acoustic modes has been studied extensively by Mack- • - 

Most of the studies on supersonic boundary layer stability were performed for flows 
over rigid surfaces. However, most aircraft structures are made up of flexible surfaces, 
therefore it is very important to study the interaction (if any) between such structu.es 
and unsteady disturbances in the boundary layer. It is known that the unsteac y pressure 
field in the boundary layer can induce significant vibrations of the flexible surface. The 
vibrating surface can radiate sound over a broad range of frecpiencies. These acoustic 
disturbances can (i) excite TS type waves in the boundary layer and (,.) poss, > y c rang 
the stability characteristics so as to enhance the instability of these waves. In turn the T 

type waves can excite vibrations of the flexible surface leading to a positive eec ac wnci 

V r Vnn Tt is therefore important to determine the coupling 
may promote an earlier transition. It is thereiore p 

between unsteady disturbances and a vibrating surface. 
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n a previous paper’ we considered the behavior of unstabie, second mode distur- 
bances m a high speed boundary layer over a flexible surface. It was shown that while the 
disturbances were unstable and exhibited substantial growth over the surface, there was 

“t‘ e e ; Cltat,0n ° ( the fleXibk "*“• This - due to the frequency range of the second 
mo e disturbances, which were too high to effectively couple with the flexible surface. 

In this paper we consider the coupling of a flexible surface with two dimensional first 
mot e i-e. TS type) disturbances. For the Mach number and Reynolds number considered 
m this paper, the two dimensional disturbance exhibits a very small growth near the 

“ ° W ’ and lhe " deCayS With * he doWnStre ™ ■ Even though the disturbance 

requency ,s an order of magnitude lower than that considered previously’, we find that 

ie two dimensional disturbance does not effectively couple with the flexible surface and 

there are no significant differences between the evolution of the disturbance over a 
flexible or rigid surface. 

However, a stronger coupling of the flexible surface with the boundary layer can be 
stained when the surface is forced by an acoustic disturbance located on the side of the 
flexible surface opposite to the boundary layer. Acoustic excitation of this sort can lead 
to large disturbances in the boundary layer. The effectiveness of the flexible surface in 
transmitting acoustic energy into the boundary layer is highly sensitive to the frequency 
of the acoustic excitation and the parameters of the flexible surface, in particular the 
‘ amping. At the present time the computational model is limited to two dimensions and 
to infinitesimal surface vibrations. However within these limitations there are indications 
that the disturbances introduced by the vibrating surface can serve to destabilize the 
flow field and that surface vibration should be considered as a potentially destabilizing 
mechanism which should be accounted for in stability studies. 

The remainder of this paper is organized as follows; first the mathematical model is 
presented, then the numerical scheme used to solve the set of partial differential equations 
is described and finally the numerical results and conclusions are presented. 
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Formulation of the model 

As shown in Fig. (1), tllr ee computation domains are being considered, the flow re- 

r >OW ““ S " rfaCe ' thS toibfe ** *»« .he no-flow region beiow the surface. 

* equations m the supersonic flow region are the two-dimensional, compress- 

" eq " ali ° ,1S - Ll * coordinate system, . and y, these equations 

can be written in conservation form as; 


Qt = F z + G„ 


where Q is the vector (p, p U , pv , E) T % p is the 
fluxes, respectively, and E is the total e: 


( 1 ) 


pu and pv are the x and y momentum 


energy per unit volume given by 


E = 2 P( u2 + v 2 ) + pc v T. 
In eq. (1), the functions F and G 


( 2 ) 


are: 


and 


pu 

F = I P u 2 +P~t xx 

P U V - T yx 

. u(E + p) ~ UT XX - VT yx - kT x 


pv 


G = I P UV - T xy 

pv*+p- Tyy I ( 3 ) 

’{ E +P)-UT yx ~VT yy ~ KTy 

where r,, are the components of the viscous stress tensor and * is the thermal conductivity, 
n ac ition to Eq. (1), an ideal gas state equation is used- 


P = pRT, 

where p is the pressure, p the density, R the 


( 4 ) 


viscosity is obtained from Sutherland’s law 


gas constant and T the temperature. The 


c, T 3 / 2 

T\Tf ’ ( 5 ) 
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with T, = 198.6 °R, and c, = 3.66x10- N sec/(.n 2 R 1 ' 2 ) tor air. 

The equation describing the motion of the flexible surface is; 

d 4 w u d 2 w dw + - 

D— + Pfh-ge + - P P 


( 6 ) 


dx 4 


where to is the plate transverse deflection, p r the mass per unit volume of the pi , 
plate thickness and 7 is the physical damping. In eq. (6), D = Mh /12(1 - ■' » 

stiffness of the plate, with M being the modulus of elasticity and „ the Poisson ratio of 
the plate material. The terms p+ and p~, are the pressure fields in the flow and no-flow 
regions, respectively. oFVom the solution of the Navier-Stokes equations we obtain P 
For p-, two cases are studied; first when a sound source is used, we consider an nnposec 

pressure field of the form, 

p" — Po d" sin(u>i t) 

Where c, and ui, are the amplitude and frequency of the acoustic source, respectively. In 
the absence of the sound source, an approximation to the solution of the wave equation 

dw 


P = Po ~ (P C )o {ft ’ 


( 8 ) 


is used This is similar to the approximation presented by Miksis and Ting» . The variables 
p- and (pc) 0 - are the free-stream pressure, density and speed of sound in the region below 
the surface. A physical situation where Eq. (8) can be used is that of flow over a use age 
surface, and Eq. (7) is a simplification of flow over a vibrating flexible surface exci ec 
by engine noise. Since the geometry used is that of a flat surface, the computations were 
limited to the case of small deflection,, (order of the plate thickness). The coupling between 
the boundary layer and the plate is achieved by imposing the vertical velocity of the p a e 
to be one of the boundary conditions for Navier-Stokes equations. Similarly, the sur ace 
motion is driven by the pressure field in the boundary layer. 

Method of solution 

The unsteady Navier-Stokes equations (eq. (1)) are solved using an explicit finite 
difference scheme. The scheme, which is a generalization of MacCormack s scheme obtaine 
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by Gottlieb and Turkel 9 , is fourth order accurate on the convective terms and second 
order accurate on the diffusive terms, and is second order accurate in time. The numerical 
scheme, applied to a one- dimensional equation of the form 


u t = F x , 


(9) 


consists of a predictor step given by 


< = K + 


At 

6A.r 


(— 7F t + 8 F I+1 -F l+2 ), 


( 10 ) 


followed by a corrector step of the form 


u 


n+1 


* I (l 

u i + u i 


At 




(ii) 


In the above equations, the subscript i denotes the spatial grid point and the superscript 
n the time level. The fourth order accuracy is obtained by alternating the scheme given 
above with its symmetric variant. In order to apply this scheme to a two-dimensional 
problem, operator splitting method is used. If L x and L y denote the solution operators 
for the one dimensional x and y problems, then the solution to eq. (1) is obtained by 


Q n+2 = L x L y L y L x Q n . 


( 12 ) 


Further details regarding the method and the advantage of fourth order schemes can be 
found in Bayliss et al. 10,11 . 

The boundary conditions employed on the surface for the Navier-Stokes equations are 

u = v = 0 and T = T w (13) 


over the rigid part of the surface, and 


u = 0, 



and T = T w 


over the flexible part of the surface. 


(14) 
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The pressure boundary conditions are as follows; over the rigid part of the surface 
the pressure is calculated using the normal momentum equation and over the flexible part 
of the surface a linear extrapolation from the interior of the domain is used to find the 
pressure. We first update the pressure to the new time level by integrating the Navier- 
Stokes equations and using dw/dt at the previous time level as a boundary condition, then 
using the new pressure field we solve the plate equation to update w. 

The inflow conditions are given by 


Q l nflou, = Qo+^[^y)e , ^ t } (15) 

where Qo is the steady state solution corresponding to a mean flow over a rigid surface, 
3? indicates the real part, (j>(y ) is the first mode disturbance eigenvector solution obtained 
from a compressible stability code, for a given set of boundary layer parameters (boundary 
layer profile and inflow Reynolds number) 12 and u ; 2 and e 2 are the frequency and amplitude 
of the inflow disturbance, respectively. The remaining inflow, outflow and upper-boundary 
conditions are the same as in Maestrello et a/. 13 . 

The plate equation is integrated using an implicit finite difference method for struc- 
tural dynamics developed by Hoff and Pahl 14 . The boundary conditions used to solve the 
plate equation are those for a clamped plate 

w = w x =0 at x = xo, x q + L. (16) 

The problem of transient oscillations is common when trying to solve the time depen- 
dent structural equations. In order to eliminate the effects of this transient on the flow 
field, the plate equation is integrated to obtain a steady state solution with the acoustic 
excitation alone. To accelerate the convergence, higher physical damping is used initially 
then it is reduced to the desired value progressively. The steady state displacement, ve- 
locity and acceleration profiles are then used as inputs to the Navier-Stokes calculations, 
which is then run to eliminate the remaining transient. The time scales of the surface 
vibration are long compared to the allowable time steps for the Navier-Stokes solver; typ- 
ically SxlO -5 sec for the structure and 10 -7 sec for Navier-Stokes. Thus there are limits 
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to the length of time for which it is currently feasible to solve the equations of the model 
We believe that the data presented here represents an approximation to the steady state 
surface response to within the accuracy of the model. 

Results and Discussion 

Numerical experiments are carried out for a supersonic laminar boundary layer with 
a freestream Mach number 2.2, a Reynolds number per meter of 5.28x10* and a total 
temperature of 311 °K. The properties of the flexible part of the surface are assumed to 
be independent of positron and are; stiffness D = 1.46 Nun, mass per umt area p,h - 2.26 
Kg/m 2 and physical damping 7 = 131.2 N-sec/m 2 . The plate is 0.264 m long and ,8. , ,m, 
thick, and is clamped between two rigid surfaces. The first seven natural frequences of the 
plate are; 45 Hz, 122 Hz, 206 Hz, 3% Hz, 501 Hz, 826 Hz and 1100 Hz. The dimensions 
o, the computational domain are 1.27 m in the downstream distance and 0.0254 m in the 
vertical distance corresponding to 20 boundary layer thicknesses. Th.s large value o t re 
vertical distance is chosen in order to calculate the freestream radiated pressure and has 
no effect on the numerical results. The number of points used ar e 301 and 201 m t u 
streamwise and vertical directions, respectively. An exponential stretching is used m the 
vertical distance in order to achieve good resolution in the boundary layer. 

The different configurations that we compute are shown m Fig. 1; a rigid surface, Fig. 
la a flexible surface clamped between two rigid ones, Fig. lb, and an acoustic source placed 
below the flexible part of the surface, Fig. 1c. The acoustic source emits plane waves at 
normal incidence to the surface at different frequencies and with a constant sound pressure 
,evel. Disturbances are introduced in the flow field at both the inflow boundary an, » 
acoustic excitation of the surface from below. A two dimensional TS type disturbance is 
introduced at the inflow with a normalized amplitude in „ of 0.08 of the freestream, and a 
normalized frequency F = 2* fu/Ul of OOxlO"* corresponding to a dimensional frequency 
of / = 4500 Hz (u .2 = 2* /). This frequency corresponds to the least stable two dunensiona 
first mode. In the definition of F, u is the kinematic viscosity and U x is the freestream 
streamwise velocity. Based on linear stability theory, this disturbance is expected to grow 
and then decay with increasing streamwise distance. In order to follow the evo ution o 
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the disturbance downstream, the mass flux disturbance level in the streamwise direction 
(pu) * calculated. First, the streamwise mass flux fluctuation in time is obtained by 


( P U )' — pu— < pu > 


(17) 


where < pu > is the computed mean obtained from integrating the data in time The 

integration is performed over one period of the sound source frequency. The disturbance 

then calculated by computing the root mean square (RMS) in time of (pu)’ and 
then integrating in y, 

RMS =< ( pu y 2 >*/2 ^ 


G(x) = J (RMS)dy. 


(19) 

Figure 2 shows the results of this calculation for a disturbance propagating over the 
rigid surface of Fig. la, and the vibrating surface of Fig. lc. The various disturbances 
used are indicated in parenthesis, where / is the TS frequency (4500 Hz). We note that 
all three computations are graphically indistinguishable upstream of the flexible surface. 
Over the rigid surface the disturbance grows and then decays as predicted from the linear 
stability theory. When a sound source emitting plane acoustic waves at a frequency of // 4 
IS used to excite the surface, the interaction between the TS disturbance and the vibrating 
surface gives rise to a small increase in the disturbance level compared to the rigid surface 
case. As the frequency of the sound source is reduced to // 8, the disturbance level clearly 
departs from that obtained over a rigid surface. The results indicate that excitation of 
sufficiently low frequency can be effective in exciting substantial vibrations of the flexible 
surface. These vibrations can significantly change the resulting evolution of the unsteady 
disturbances in the flow field. In the absence of the sound source, the disturbance level 
over the flexible surface is the same a, that over a rigid surface indicating a weak coupling 
between the surface and the disturbance. 

The effect of sound source frequency was investigated in the absence of the TS type 
waves and the results are shown in Fig. 3. For the same sound pressure level of the 
acoustic excitation Fig. 3 shows that the the disturbance level introduced by the flexible 
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surface increases as the frequency decreases. For high frequency excitation the disturbance 
level in the boundary layer is nearly uniform over the flexible surface. As the frequency 
is reduced the disturbance level grows and is concentrated near the leading edge of the 
flexible surface. However, in all the cases the disturbance level decays over the trailing 
edge and downstream of the flexible surface. 

In order to further explain the differences shown in Fig. 2 between the // 4 and //8 
frequencies, a comparison of the disturbance level obtained for a rigid surface, a flexible 
surface excited by sound alone and a flexible surface excited by a combination of a TS type 
waves and sound, was made. For a sound source frequency of //4, the disturbance level due 
to the sound source alone is everywhere less than that of the TS type disturbance over the 
rigid surface. Therefore, the level of the combined disturbance over the flexible surface is 
only slightly higher than that of a rigid surface as is shown in Fig. 4a. However, when the 
frequency of the sound source is reduced to //8 the level of the disturbance generated by 
the vibrating surface is higher than that of the TS type disturbance over the rigid surface. 
This results in a higher level of the combined disturbances over the flexible surface, Fig. 
4b. It is interesting to note that near the leading edge of the flexible surface, the level of 
the combined disturbances is higher than that of the sound source alone; however, near the 
trailing edge the two curves coincide. This is an indication that the two disturbances do 
not simply add up, but that there are regions of constructive and destructive interference 
between the two waves. 

Figure 5 shows the instantaneous pressure distribution in the flow field over the rigid 
and flexible surfaces. The pressure distribution over the rigid surface is generated by 
the instability waves, Fig. 5a. However, when the same disturbance propagates over a 
vibrating flexible surface it interacts with the pressure disturbance radiated by the surface. 
Outside the boundary layer, the radiated pressure curves in the direction of the flow and 
propagates along the Mach line (at a Mach angle of 27° ). Figure 5b shows the interaction 
between these two waves when the frequency of the acoustic source is // 4. This excitation 
frequency (1125 Hz) is near the seventh natural frequency of the surface, and therefore 
its response and near field radiation patterns are dominated by the seventh mode. This 
seventh mode patterns are shown in the pressure distribution of Fig. 5b. Decreasing the 
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sound source frequency to // 8 (563 Hz) leads to the excitation of lower modes on the 
surface. The pressure field shown by Fig. 5c corresponds to that radiated by a third 
mode. The increase in vibration level and the mode shape of the flexible surface result 
in greater coupling between flow and structure which explains the increase in disturbance 
level for the // 8 case. One should note that Figs. 5b-c show only a small fraction of 
a wavelength of the acoustic wave radiated by the surface because of the small vertical 
dimension (0.0254 m) compared to the acoustic wavelength (0.61 m for // 8). 

Figure 6 shows a comparison of the instantaneous velocity profiles at the center of the 
flexible and rigid surfaces for the three cases of Fig. 2. When the frequency of the acoustic 
source is f j 8, the velocity profile has an inflection point which suggests a more unstable 
flow field (inflectional instability) as compared to the other two profiles. The variation of 
the velocity profile with time for the case when the frequency of the acoustic source is / /S, 
is shown in Fig. 7. As the surface motion goes through one cycle (based on // 8, period 
T), the velocity profile varies between the upper and lower profiles. The change in the 
velocity profile leads to a change in boundary layer thickness. A thicker boundary layer 
is associated with the upper, less stable profile while a thinner boundary layer is obtained 
for the lower profile. This result indicates that the flexible surface acts like a piston, with 
blowing and suction phases. 

The vorticity contours of the flow field generated by the TS disturbance over the 
rigid surface are shown on Fig. 8a, and those generated by the interaction of the TS 
disturbances with a vibrating flexible surface excited by a sound source of frequency // 8 
are shown on Figs. 8b-c. As the flexible surface moves upward, the vorticity contours 
cluster near the leading edge of the flexible surface indicating an increase in disturbance 
level in that region, Fig. 8b, which is in agreement with Figs. 3 and 4. In addition, the 
vorticity is lifted away from the wall, leading to an increase in boundary layer thickness. 
Near the leading edge of the flexible surface, there are indications of possible vortex roll 
up which can be a precursor to transition. When the surface is moving downward, the 
vorticity contours are concentrated near the flexible surface, Fig. 8c. The boundary layer 
becomes thinner and there is no indication of vortex roll up. 


10 



Conclusions 


Based on the results presented here the following conclusions can be made: 

i) The coupling between TS type waves and a flexible surface is not significant. This 
result is caused by the fact that the frequencies of the TS type waves are too high 
to effectively excite the surface. 

ii) Acoustic excitation of the flexible surface can effectively excite vibrations of the 
surface, provided the frequency is sufficiently low (for the given set of structural 
parameters). In this case the resulting surface vibration can significantly enhance 
the level of unsteady disturbances in the flow and change the nature of the flow 
field. In particular inflection points can develop in the velocity profile and the 
incipient formation of vortex rolls can be seen. 

iii) The effect of the surface vibrations can be transmitted outside of the boundary 
layer as acoustic radiation. 
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Fig. 1: Computational Domains, (a) rigid surface, (b) flexible surface without 
sound, (c) flexible surface with sound. 
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Fig. 2: Comparison of streamwise mass flux disturbance level generated by the 
interaction of a two dimensional disturbance with a rigid surface and a vibrating 
flexible surface at different downstream locations. 
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Fig. 3: Effect of sound source frequency on the streamwise mass flux disturbance 
level at different downstream locations. 
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Fig. 4: Comparison of disturbance levels generated by a TS type waves over a 
rigid surface to that generated by a vibrating flexible surface (no TS waves) and 
to that generated by the interaction of the TS waves with the vibrating suiface. 
(a) Sound source frequency of // 4, (b) Sound source frequency of // 8. 
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>: Instantaneous pressure distribution in the flow-field, (a) TS type waves 
rigid surface, (Id) TS type waves over a vibrating flexible surface, sound 
frequency // 4, (c) sound source frequency f/8. 
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Fig. 6: Comparison of instantaneous velocity profiles at the center of a rigid 
surface and a vibrating flexible surface. 
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Fig. 7: Time variation of the velocity profile at the center of the vibrating 
flexible surface, sound source frequency // 8 (T is the oscillation period). 
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Fig. 8: Instantaneous vorticity contours in the flow field generated by TS waves, 
(a) Over a rigid surface. A vibrating flexible surface with sound source frequency 
// 8 (b) surface moving upward (c) surface moving downward. 
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